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ABSTRACT 

Many scenarios for the origin of the chemical anomalies observed in globular clusters 
(GCs; i.e., multiple populations) require that GCs were much more massive at birth, 
up to 10 — 100 X, than they are presently. This is invoked in order to have enough 
material processed through first generation stars in order to form the observed numbers 
of enriched stars (inferred to be second generation stars in these models). If such 
mass loss was due to tidal stripping, gas expulsion, or tidal interaction with the birth 
environment, there should be clear correlations between the fraction of enriched stars 
and other cluster properties, whereas the observations show a remarkably uniform 
enriched fraction of 0.68 ± 0.07 (from 33 observed GGs). If interpreted in the heavy 
mass loss paradigm, this means that all GGs lost the same fraction of their initial 
mass (between 95 — 98%), regardless of their mass, metallicity, location at birth or 
subsequent migration, or epoch of formation. This is incompatible with predictions, 
hence we suggest that GGs were not signihcantly more massive at birth, and that the 
fraction of enriched to primordial stars observed in clusters today likely reflects their 
initial value. If true, this would rule out self-enrichment through nucleosynthesis as a 
viable solution to the multiple population phenomenon. 

Key words: galaxies - star clusters. Galaxy - globular clusters 


1 INTRODUCTION 

Most scenarios put forward to explain the observed abun¬ 
dance anomalies in globular clusters (GCs) invoke some form 
of self-enrichment, where stars produce “enriched material” 
through nucleosynthesis (i.e. material that displays the ob¬ 
served abundance variations in light elements) that pollutes 
other stars within the same cluster. In the popular “AGB 
scenario” (e.g., D’Ercole et al. 2008), AGB stars, formed 
in a first generation, shed enriched material which, when 
mixed with primordial material, form a second generation 
(or more) of stars that show the observed chemical anoma¬ 
lies. Due to AGB star lifetimes, this is expected to operate 
over 30 — 200 Myr timescales. The “Fast Rotating Massive 
Star” scenario (FRMS - e.g. Decressin et al. 2007; Krause 
et al. 2013) also invokes multiple star-forming epochs, al¬ 
though the timescales are significantly shorter. This scenario 
uses the ejecta of rapidly rotating massive stars, when mixed 
with primordial material left after the formation of the first 
generation, to form subsequent generations. The timescales 
associated with the FRMS scenario are < 10 Myr. 

In these scenarios, stars with primordial abundances 
(i.e. similar to that shown in halo field stars) are interpreted 
to be from the first generation (FG) while stars that are en¬ 


riched in some light elements (i.e., Na, Al, He) and depleted 
in others (C, O) are interpreted to be from a second genera¬ 
tion (SG). We note that the observations only show whether 
a star is chemically depleted/enriched or not, whereas the 
labels of first or second generation are interpretations, as 
there is no direct evidence for the majority of GCs of mul¬ 
tiple star-forming events (see the discussion in Chantereau 
et al. 2015). 

One problem that the above scenarios have is that ob¬ 
servations show that there are approximately equal numbers 
of primordial (assumed to be FG in the above scenarios) and 
enriched (assumed to be SG stars) stars within GCs (e.g., 
Carretta et al. 2009). Since the polluting stars only make up 
a very small fraction of the total initial mass of the cluster 
(3 — 8% - e.g., de Mink et al. 2009) the current observed pri¬ 
mordial stars (and their associated higher mass stars that 
have already burned out) are not enough to have produced 
the required amount of enriched material. This is often re¬ 
ferred to as the “mass budget problem”. 

A commonly invoked solution to this fundamental prob¬ 
lem of self-enrichment scenarios is to have the clusters lose 
large fractions (> 90%) of their FG stars due to tidal field 
of the Galaxy (e.g., D’Ercole et al. 2008; Conroy 2012) or 
gas expulsion after the SG has formed (after ~ 10 Myr for 
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the FRMS scenario or ~ 100 Myr for the AGB scenario). 
This itself would need some fine tuning to result in all clus¬ 
ters having similar enriched stellar fractions, fenriched, from 
an initial fenriched^S 0-05 to the current fenriched= 0.5 — 0.8, 
implying that every cluster lost 95 — 99% of their initial 
FG stars. While finely tuned, this is, at least in principle, 
feasible. 

Observations of the ratio of field and globular cluster 
stars in dwarf galaxies, however, have shown that GGs in 
those systems could not have been more than a factor of 

5 — 10 (and potentially substantially less) more massive at 
birth than they currently are (Larsen et al. 2012; 2014a), 
calling into question the proposed scenarios. Schiavon et 
al. (in prep.) have found, based on the recent APOGEE 
survey, that GGs in the bulge of the Galaxy could only 
have had initially ~ 8 times more primordial stars than 
enriched stars (and likely much less), again in contradiction 
with the model requirements. The Larsen et al. and Schiavon 
et al. constraints are quite complimentary, as they sample 
very different environments (dwarf galaxies and the bulge 
of the Milky Way, as well as different metallicity regimes, 
[Fe/H]< —2 and —2 < [Fe/H] < —1, respectively). 

There are other independent tests that can be carried 
out on the GG population of the Milky Way in order to 
test whether the requirements of the self-enrichment models 
can be satisfied. One way is to look for the effects of tidal 
stripping on the multiple populations within current GGs. 

If mass loss is due to tidal stripping] from the host 
galaxies, the dissolution time of a cluster is linearly propor¬ 
tional to the Galactocentric distance (assuming a flat rota¬ 
tion curve - e.g., Baumgardt & Makino 2003; Kruijssen & 
Mieske 2009; Lamers et al. 2010). Hence, for a given mass, 
clusters in the inner regions of the Milky Way should lose 
their stars more rapidly than clusters in the outer regions. 
Because the scenarios invoke the preferential loss of FG 
stars, fenriched would increase as the cluster dissolves. The 
result should be that clusters in the inner parts of galaxies 
should have a higher fraction of SG stars, while GGs in the 
outer parts of the Galaxy should be dominated by their FG 
stars. The same applies to GGs formed within dwarf galax¬ 
ies, as their tidal fields are substantially weaker than the 
inner Galaxy. 

This can be tested by looking for a trend of fenriched with 
Galactocentric distance. Di Criscienzo et al. (2011) studied 
the Galactic cluster, NGG 2419, a massive (~ 10® Mq) GG 
that is one of the most distant (-^ 90 kpc) of the Milky Way 
population. In the strong cluster dissolution paradigm, in¬ 
voked in the self-enrichment scenarios discussed above, due 
to its high mass and large Galactocentric distance, this clus¬ 
ter is not expected to have lost many FG stars, so its current 
ratio should be close to the initial value. Contrary to ex¬ 
pectations, Di Criscienzo et al. (2011) found that enriched 
stars make up > 30% of the population, which has been 
spectroscopically confirmed (Mucciarelli et al. 2012; Cohen 

6 Kirby 2012). The authors recognised this discrepancy, and 
as a potential solution, suggested that the orbit of NGG 2419 


^ We use the term “tidal stripping” in a general way to include 
the loss of GC stars due to two-body relaxation in a tidal field 
(where stars escape more easily due to the tidal boundary) as well 
as due to the expansion of the cluster over the tidal boundary. 


was extremely elliptical, with a peri-centre passage of only 
~ 11 kpc. However, due to its high mass even at this dis¬ 
tance, NGC 2419 would not be expected to lose a large frac¬ 
tion of its stars (e.g., Baumgardt & Makino 2003). Hence 
the high enriched fraction of stars within this cluster re¬ 
mains unexplained. 

We can also search for relations between fenriched and 
metallicity, as it is low-metallicity clusters (in the outer 
Milky Way) that are expected to be accreted from dwarf 
galaxies. This will test whether the birth environment of 
GCs is responsible for the required mass-loss to make self¬ 
enrichment scenarios viable, although we note that current 
estimates do not find the heavy mass loss required for the 
majority of GGs (e.g., Kruijssen 2014; 2015). 

While any individual GG may have an anomalous orbit 
(i.e. highly elliptical), it would require extreme fine tuning if 
other distant GGs were found that had significant fractions 
of enriched stars. This would require the GG population has 
a whole to have a highly radially anisotropic velocity dis¬ 
tribution, which is not consistent with observations (see the 
discussion in Vesperini et al. 2003). In the present paper we 
use data from the literature on the fraction of enriched stars 
within particular GCs based on large spectroscopic datasets, 
and supplement this with studies that use HST photome¬ 
try in the appropriate filters where the different populations 
can be distinguished. We then search for the expected radial 
trend, cluster mass relation, or metallicity influence, which 
are required in self-enrichment scenarios. The expectation 
is that outer halo GCs, higher mass GCs, or low metallicty 
GCs should have a drastically reduced fraction of enriched 
stars (reflecting their initial fractions). 

Finally, we can search for a relation between fenriched 
and cluster mass. In the one-dimensional hydrodynamical 
simulations performed by D’Ercole et al. (2008) and Ves¬ 
perini et al. (2010), the initial fraction of fenriched is a strong 
function of mass, with lower mass (10® Mq) clusters initially 
having fenriched= 0.01 and higher mass clusters (-^ 10® Mq) 
initially having fenriched> 0.07 — 0.13, depending on the ini¬ 
tial radius. Hence, we may expect that some remnant of this 
relation has been maintained during the subsequent evolu¬ 
tion of the GCs. In the present work we adopt an initial 
value of /enHched = 0.05 from the models of e.g., Vesperini 
et al. (2010). 

This paper is structured as follows. In §[5] we introduce 
the datasets used in the present work while in §[3]we investi¬ 
gate the expected relations between fenriched and other clus¬ 
ter properties if the clusters have undergone large amounts 
of mass loss during their evolution. In §|4]we discuss our re¬ 
sults in light of the self-enrichment scenarios as well as how 
our results relate to others in the literature. Finally, in § [S] 
we present our conclusions. 


2 OBSERVATIONS 

The primary dataset used in the present work is taken from 
Carretta et al. (2009; 2010a). The authors used a large spec¬ 
troscopic database of red giant branch stars (RGBs) in a 
sample of 19 Galactic GCs, and determined the fraction 
of primordial (P), intermediate (I) and extreme (E) stars., 
where I and E stars show enrichment patterns not observed 
in halo held stars. For the present work, we sum the frac- 
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tions of the I and E populations, as we are only concerned 
about the fraction of enriched stars, regardless of the level 
of enrichment that they display. 

To this dataset we include results from a number of re¬ 
cent photometric and spectroscopic samples. We add the 
following spectroscopic datasets: NGC 362 (Carretta et 
al. 2013); NGC 1851 (Carretta et al. 2011); NGC 4833 
(Carretta et al. 2014a); NGC 5286 (Marino et al. 2015); 
NGC 6093 (Carretta et al. 2015), NGC 6121 (Carretta et 
al. 2010a); NGC 6752 (Carretta et al. 2012); and NGC 6864 
(Kacharov et al. 2013). We have also included the spec¬ 
troscopic results of NGC 5272 and NGC 6205 (Sneden et 
al. 2004; Cohen & Melendez 2005) which have been anal¬ 
ysed in terms of their P, I and E assignment by Carretta et 
al. (2009). 

We include the spectroscopic results on the distant and 
massive GC, NGC 2419, from Cohen & Kirby (2012). If we 
use [Mg/Fe] to split enriched vs. primordial stars we ob¬ 
tain fenriched = 0.61. If, instead, we use [Na/Fe] as the dis¬ 
criminator we obtain fenriched = 0.84. For the present work 
we simply take the mean of these two values and adopt 
fenriched = 0.73. We uote that in general [Mg/Fe] is not 
a good tracer of chemical anomalies, as [Mg/Fe] does not 
vary in many clusters (c.f. Bastian et al. 2013b) that are ob¬ 
served to host multiple populations, hence [Na/Fe] is likely 
more representative for the cluster. This is consistent with 
the lower limits inferred spectroscopically by Mucciarelli et 
al. (2012) and photometrically by di Criscienzo et al. (2011). 
These are lower limits, as the two studies were most sensi¬ 
tive to extremely enriched (either depleted in Mg or enriched 
in He) stars, and likely missed the intermediate population 
stars. 

Similarly, the results of NGC 6266 (M62) from Yong et 
al. (2014) are included, although this point has relatively 
large errors as it is based on only 7 stars (4 enriched, 3 
primordial, resulting in fenriched = 0.57). Just taking Pois- 
sonian errors leads to an uncertainty of ~ 0.3. This result 
agrees with the photometric study of Milone (2015), who 
found a lower limit of fenriched> 0.2. This result was a lower 
limit as the author detected two clearly distinguishable se¬ 
quences along the main sequence of the cluster, however, 
the main sequence associated with the “primordial” popu¬ 
lation was broader than the photometric errors (consistent 
with a spread in helium). Hence, the stars associated with 
this main sequence likely represent a mixture of primordial 
and intermediate enriched stars. Additionally, the photom¬ 
etry presented in Milone (2015) was more sensitive to He 
spreads than spreads in light element abundances and some 
clusters with relatively large Na and O spreads display small 
He spreads (e.g., Bastian et al. 2015). 

The use of high precision photometry of GC stars with 
specific filters that are sensitive to light element abundance 
variations has the potential to significantly refine the esti¬ 
mates of fenriched based OH spectroscopic measurements (e.g., 
Milone et al. 2015). 

We note that we have not homogenised the data (i.e. 
applied consistent NLTE corrections to the Na abundances), 
although the majority of the datasets come from the same 
team, which use consistent techniques. 

We also use the results based on HST photometry of 
NGC 6362 (Dalessandro et al. 2014) and NGC 7089 (Milone 
et al. 2015) 


Finally, we include the estimate of D’Antona & 
Caloi (2008) for NGC 6229. This result is based on the Hor¬ 
izontal Branch morphology from HST photometry. Since it 
relies on detailed stellar evolutionary modelling, instead of 
coming directly from the observations, we consider this to 
be a relatively uncertain data point. However, it follows the 
overall trends for the other clusters in our sample. 

In several cases, photometric and spectroscopic datasets 
exist for the same cluster, and estimates of fenriched have 
been made using both datasets. Overall, the agreement be¬ 
tween the techniques is good, i.e., generally within the re¬ 
spective error bars (e.g., NGC 6121 - Carretta et al. 2010a 
and Milone et al. 2014 or NGC 2808 - Carretta et al. 2010a 
and Piotto et al. 2007). However, we note that as photo¬ 
metric accuracy improves along with imaging in filters that 
are good tracers of chemical anomalies, the picture is be¬ 
coming a bit more complicated. For example, the red main- 
sequence identified in NGC 2808 by Piotto et al. (2007) 
which is thought to correspond to the “primordial” stars 
identified spectroscopically, actually consists of three sub¬ 
populations (Milone et al. 2015). Hence, it is unclear which 
of the three sub-populations of the “red-MS” actually cor¬ 
responds to the primoridal population. This has also been 
noted a study of NGC 6752 by Milone et al. (2013), who 
found three populations based on photometry, which did 
not fit neatly into the P, I and E spectroscopic framework 
of (Carretta et al. 2009a). However, even here, there was a 
reasonable agreement in distinguishing between P and I(-l-E) 
stars using photometric and spectroscopic methods. Future 
studies, using a homogenised analysis of a full sample of GC 
with the appropriate filter combination (e.g., with the HST 
UV Legacy Survey of GCs - Piotto et al. 2015) will compli¬ 
ment and extend the current study. 

Based on the comparison between spectroscopic and 
photometric divisions of the sub-populations, some of the 
estimated fenriched from spectroscopy may be lower limits. 

We have intentionally excluded GCs that are commonly 
referred to as members of the Sagittarius Dwarf Galaxy as 
their current position within the Galaxy is not represen¬ 
tative of the tidal field experienced by the clusters over a 
majority of their lifetimes. We note that this removes M54 
and NGC 4590; both which show large enriched populations 
(Carretta et al. 2010b), but also the low mass clusters Terzan 
7, 8 and Palomar 12 which hosts small or non-existent en¬ 
riched populations according to spectroscopic studies (e.g., 
Carretta et al. 2014b). We note that Terzan 7 and Palomar 
12 are much younger than typical GCs (< 9 Gyr - Teaman 
et al. 2013) hence age may play a role in whether chemical 
anomalies are present in clusters, although Terzan 8 appears 
to be similar in age to the bulk of Galactic GCs (Marm- 
Franch et al. 2009). Future imaging of these clusters should 
provide a better estimate of fenriched for these clusters. 

Additionally, we have excluded GCs with evidence of 
significant spreads in [Fe/H], such as w-cen, M54, M22 and 
M19. This is because clusters with an intrinsic iron spread 
also display an s-process element bimodality, with the s-rich 
population also enriched in iron. Each s-process groups has 
its own C-N and Na-0 anti-correlation (s-rich stars are also 
on average richer in C, N, and Na with respect to s-poor 
stars). In theses cases it is not trivial to distinguish between 
P, I, and E stars, since the two Na-0 anti-correlations over¬ 
lap (e.g., Marino et al. 2011). 
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This results in a sample of 33 GCs with full information 
available. 

We show the resulting trend between the fraction of en¬ 
riched stars, fenriched, and Galactocentric distance, GG mass 
and metallicity in Fig. [T] Gluster masses were taken from 
Kruijssen & Mieske (2009) or Gnedin & Ostriker (1997), 
while Galactocentric distances and [Fe/H] were taken from 
Harris (1996, 2010). Throughout this work the fractions re¬ 
ferred to are that of the long lived low mass stars, in order 
to be able to compare with observations directly. 

The data show a remarkably constant enriched frac¬ 
tion, with little scatter, fenriched = 0.68 ± 0.07 (the mean 
and standard deviation of the sample). This finding is simi¬ 
lar to that found in Carretta et al. (2009 - their Fig. 10), but 
now with a larger GG sample. Garretta et al. (2009) have 
estimated the expected scatter, based on the observed num¬ 
bers of stars and errors associated with assigning stars to 
the P, I or E sub-populations, of (T(fenriched) — 0.07 for their 
sample (making up 19 out of 33 GCs in our sample). Hence, 
the observed scatter can be nearly entirely explained within 
the expected error budget, leaving little room for correla¬ 
tions with other variables. Least squares fit show that the 
fenriched is cousisteut with uo relation with Galactocentric 
radius, mass or metallicity of the clusters. 

From the data directly we note that the similar value of 
/enriched obscrved in all elusters argues against any environ¬ 
mentally dependent dynamical process to signifieantly alter 
the initial value, as a mueh larger scatter would be expeeted. 
If interpreted as being due to mass loss, the data suggest that 
all elusters lost nearly identical fractions of their initial mass 
(between 95 — 99%) regardless of their mass, metallicity, or 
location at birth or subseguent migration. 

3 EXPECTATIONS OF CLUSTER 

DISSOLUTION 

3.1 Evolution of Clusters in a Tidal Field 

The dissolution of globular clusters due to the combined ef¬ 
fects of stellar evolution, two body relaxation within a tidal 
field and external shocks (such as disk/bulge shocks or close 
GMG passages) has been extensively studied, both numer¬ 
ically (e.g., Baumgardt & Makino 2003) and analytically 
(e.g., Lamers et al. 2010 and references therein). A basic 
result of these investigations is that for a given mass and 
density, clusters in stronger tidal fields dissolve more rapidly 
than clusters in weaker tidal fields. 

The mass budget problem within self-enrichment sce¬ 
narios has been recognised from an early stage, hence mod¬ 
els, such as D’Ercole et al. (2008), have long invoked cluster 
dissolution as a potential solution. These models have max¬ 
imised the role of dissolution by adopting extreme param¬ 
eters. Namely, they have assumed strong tidal fields (e.g., 
at 4 kpc from the Galactic centre), very diffuse clusters 
(half mass radii of 25 pc - in contrast with young massive 
clusters and the present properties of GGs - e.g., Portegies 
Zwart et al. 2010), that the clusters are tidally limited from 
birth (so that any expansion leads to FG star mass loss), 
that the clusters were initially mass segregated (to max¬ 
imise cluster expansion) and that the SG stars were much 
more concentrated initially than the FG stars (see Khalaj & 


Baumgardt 2015 for a further discussion of these points). If 
these extreme assumptions are relaxed, the models are not 
able to explain the observed ratios of enriched and primor¬ 
dial stars (D’Ercole et al. 2008). 

Baumgardt & Makino (2003) have carried out a large 
suite of N-body simulations, which include stellar evolu¬ 
tion (and the resultant stellar evolution driven expansion), 
of clusters of various masses and at different Galacocentric 
distances. The authors find that the dissolution time of a 
cluster is directly proportional to the Galactocentric radius 
(tdis oc Rgc), for a given mass and orbital eccentricity. This 
result is also expected from analytical theory of cluster dis¬ 
solution (e.g., Kruijssen & Mieske 2009; Lamers et al. 2010, 
and references therein). They also found that higher mass 
clusters have longer dissolution timescales, with tdia oc 

Since FG stars are generally thought to form in a more 
extended distribution than SG stars (e.g., Lardo et al. 2011; 
although see Larsen et al. 2015 for a counter-example), the 
loss of stars due to cluster dissolution will preferentially 
remove FG stars. Therefore, clusters that have undergone 
more mass loss (i.e., the loss of stars) should have a higher 
fraction of enriched stars, fenriched. Since the dissolution 
timescale (which is proportional to the rate of mass loss 
- e.g., Lamers et al. 2010) is a linear function of Galactocen¬ 
tric distance, this naturally suggests that fenriched should be 
a linear function of Galactocentric distance. 

Most of the GGs that have had their abundances anal¬ 
ysed to date are relatively high mass clusters (> 2x 10® Mq), 
where the tidal stripping is not expected to remove large 
fractions of the cluster mass, even at relatively small Galac¬ 
tocentric distances (-^ 5 — 10 kpc - e.g., Kruijssen 2014). 
Hence, most of the GGs studied spectroscopically to date 
have dissolutions timescales well in excess of a Hubble time 
(e.g., Baumgardt & Makino 2003). We note that in the 
D’Ercole et al. (2008) model, much of the mass loss is due 
to stellar evolutionary expansion and not two-body relax¬ 
ation. This process (stellar evolutionary driven expansion) 
is included in the Baumgardt & Makino (2003) simulations, 
but is expected to be weaker due to the less extreme initial 
conditions (half-light radii, mass segregation) assumed by 
Baumgardt & Makino (2003). 

Even in the case of mass-loss driven by stellar evolu¬ 
tionary expansion, the amount of mass lost is expected to 
be inversely proportional to GG mass and Galactocentric 
radius, although it may be somewhat different in magni¬ 
tude relative to two-body relaxation driven tidal stripping. 
The effects of primordial mass segregation, and its role in 
cluster dissolution through stellar evolutionary induced clus¬ 
ter expansion, has been studied numerically by Haghi et 
al. (2014). The authors carried out a suite of N-body simu¬ 
lations of clusters in a Milky Way type potential, for clusters 
with and without mass segregation. As expected, they find 
that mass-segregated clusters dissolve more rapidly (by a 
factor of ~ 4), due to stellar evolutionary driven expansion. 
However, they also find that the dissolution timescale for 
mass segregated and non-mass segregated clusters depends 
on Galactocentric radius in the same way, with Tdis oc Rccy 
and a = 1.1 — 1.3. This is similar to that found by Baum¬ 
gardt & Makino (2003). Hence, whether cluster mass-loss is 
driven mainly by two-body relaxation or stellar evolutionary 
driven expansion in a tidal field, the fraction of mass loss 
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Figure 1. The fraction of enriched stars as a function Galactocentric distance (left panel), mass (centre panel), and [Fe/H] (right panel). 
Blue (circles), green (triangles), and red (squares) symbols represent GCs with current masses < 3 X IO^Mq, 3 — 10 x IO^Mq, and 
> 10 ®Mq, respectively. The horizontal dashed lines represent the initial value in the AGB or FRMS scenario, where the SG stars only 
make up a small fraction of the total mass of the cluster, and subsequent fenrichedi if 90> 95 or 98% of the FG mass has been removed 
(from bottom to top, respectively). The solid line in the left panel shows the expected relation if strong mass loss due to tidal stripping 
acted on the clusters, preferentially removing only FG stars. This has been scaled to match the observations at 3 kpc, whereas the actual 
cluster mass loss is expected to be much weaker. The same relation is expected if stellar evolutionary driven expansion is the cause of the 
mass-loss through mass segregation. The solid line in the centre panel is the prediction of Khalaj &: Baumgardt (2015) if the mass loss 
was due to gas expulsion. A representative uncertainty in fenriched is shown in the left panel. The grey shading denotes the mean and 
standard deviation of the observed GCs. The Spearman Rank correlation coefficients are given in each panel, no significant correlation 
between fenriched and Galactocentric distance, mass or metalliicty is found. 


(and subsequently fenriched for the models discussed here) 
should be a linear function of Galactocentric distance. 

Although studies have shown that GCs are not ex¬ 
pected to have lost large fractions of their initial mass due 
to tidal stripping (e.g., Baumgardt & Makino, Kruijssen & 
Mieske 2009; Lamers et al. 2010; Kruijssen 2014), we are 
not interested in testing the absolute timescales for cluster 
dissolution in the present work. Rather we will focus on its 
expected dependence with Galactocentric distance. As such, 
we assume that cluster dissolution operates in such a way 
as to remove 98% of the FG stars at a radius of 3 kpc, and 
hence match the mean of the observations at that radius. In 
Fig.m left panel, we show the expected decline in fenriched as 
a function of Galactocentric distance, where the fraction of 
FG stars lost depends linearly on distance. This technique 
also removes much of the dependence of GG dissolution on 
cluster mass, unless there is a strong relation between mass 
and Galactocentric data. In Fig.[T]the GCs are colour-coded 
in three mass bins (see caption for details), and no strong 
dependence on cluster mass is apparent. 

Such tidal stripping should also depend on the mass of 
the GC, with higher mass GCs undergoing less mass loss 
(tdis oc - e.g., Baumgardt & Makino 2003; Lamers et 
al. 2010). Hence, we would expect that higher mass GCs 
should have low fenriched values. In the centre panel of Fig.[T] 
we show the observed relation between fenrichedand cluster 
mass. No correlation is identified, in contrast with model 
predictions. 

Additionally, we note that if clusters were substantially 
larger in the past, the amount of mass lost through two-body 
relaxation and tidal stripping would be much less, and even 
tidal shocks (bulge, disk or GMCs) would not be expected 
to remove much mass. 


We have explicitly assumed that only FG stars are lost 
(which maximises the efficiency of tidal stripping to bring 
the initial fenriched up to the current values), although even 
in the outer regions of GCs today, enriched stars are only 
slightly more centrally concentrated than stars with pri¬ 
mordial abundances (e.g., Lardo et al. 2011; Vanderbeke et 
al. 2015), hence in practice we would expect some fraction 
of SG stars to be lost in addition to FG stars, in the self¬ 
enrichment plus mass-loss models. 

In the above estimates we have assumed (as in previ¬ 
ous works, e.g., D’Ercole et al. 2008) that a cluster’s current 
orbit is representative of the average tidal field experienced 
by the GC over its lifetime. However, this may not be valid 
as many, especially the low metallicity, GCs are thought 
to have been accreted from dwarf galaxies (e.g., Brodie & 
Strader 2006). As such, we may expect a relation between 
the metallicity of a GC (as a proxy of whether it formed in- 
situ or was later accreted from a dwarf galaxy) and fenriched. 
The tidal field of a dwarf galaxy is substantially weaker than 
that of the inner Milky Way, so whether a GC spent most of 
its life orbiting at large Galactocentric radii, or was accreted 
from a low-mass dwarf (which in turn corresponds, on av¬ 
erage, to the GCs at large Galactocentric radii), we would 
expect these clusters to have lost less of their FG stars than 
inner Galaxy GCs. 

In the right panel of Fig. [T] we show fenriched against 
[Fe/H] for our cluster sample. As was found for Galactocen¬ 
tric distance and GC mass, there is no clear relation between 
the fraction of enriched stars and metallicity, so the birth lo¬ 
cation or subsequent migration has not influenced fenriched- 

Hence, we conclude that if mass-loss is driving the ob¬ 
served value o//enriched, then the mass loss is independent 
of Galactocentric radius, metallicity, current mass and birth 
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environment (and subsequent migration), in contradiction 
with all mass loss mechanisms proposed so far. 

3.2 Gas Expulsion 

Khalaj & Baumgardt (2015) have investigated, through N- 
body simulations, whether the observed high fractions of 
enriched stars could be due to mass loss induced by gas 
expulsion (following on the work of Decressin et al. 2010). 
If the SG stars form from a gas reservoir that contains a 
substantial fraction of the total mass of the cluster (i.e. sim¬ 
ilar to the total mass of FG stars), and if this gas is removed 
rapidly from the cluster, the resultant expansion of the clus¬ 
ter can result in the loss of large fractions of FG stars (we 
again are assuming to be less centrally concentrated than 
SG stars). 

Note that this is not the gas expulsion after the forma¬ 
tion of the FG, as that would remove the stars too early, 
and not allow them to contribute to the material needed to 
form the SG stars. 

Khalaj & Baumgardt (2015) find that in order to re¬ 
produce the observed fenriched, a considerable amount of fine 
tuning is necessary. Only a small part of parameter space 
was able to lead to clusters that 1) survived and 2) had 
fenrichedbetween 0.5 and 0.8. However, for the surviving clus¬ 
ters, their models would predict a strong correlation between 
the mass of the cluster and fenriched. This is due to the fact 
that in order to have large fenriched ratios, GCs need to lose 
a large fraction of the total mass, resulting in lower mass 
clusters on average. 

Their prediction (approximate linear ht to the models 
presented in their Fig. 6) is shown compared to the data 
in the middle panel of Fig. [T] We see that the expected re¬ 
lation between cluster mass and fenriched is not followed by 
the observations. Hence, we conclude, like Khalaj & Baum¬ 
gardt (2015) that gas expulsion following the formation of 
the SG is unlikely to be a viable mechanism to remove large 
fractions of FG stars. 

A similar argument could be made for the “cruel cra¬ 
dle” mass loss mechanism (c.f., Kruijssen 2014), where the 
gas rich environment (made up of GMCs) where clusters 
form, also rapidly destroys clusters (due to GMG interac¬ 
tions). This mechanism is expected to operate too quickly 
for the AGB scenario, as it is most pronounced in the hrst 
10-20 Myr of a cluster’s life. Regardless, clusters with high 
fenriched would be expected to be lower-mass, on average, as 
they would be the clusters that have undergone the most 
mass loss. Hence, if this mechanism was responsible, we 
would expect strong trends of fenriched with cluster mass, 
and potentially metallicity as this may reflect the birth en¬ 
vironment, neither of which are seen in the observations. 


4 DISCUSSION 

4.1 Constraining the Mass Loss from GCs 

The observations show that fenriched is remarkably constant 
(0.68 ± 0.07). If the current value of fenriched is interpreted 
as being due to mass loss, the data suggest that all clusters 
lost nearly identical fractions of their initial mass (between 
95 — 99%) regardless of their mass, metallicity, or location 


at birth or subsequent migration, requiring a high level of 
fine-tuning. 

We have searched for dependencies of fenriched on other 
cluster properties, which are expected to be present if clus¬ 
ters have lost large fractions of their initial masses. We find 
that it is independent of Galactocentric distance, and does 
not follow the expected trend of decreasing fenriched with 
Galactocentric radius that would indicate that a significant 
amount of mass loss has occurred due to tidal stripping 
(either through two-body relaxation or cluster expansion 
driven by stellar evolutionary mass loss in a mass segre¬ 
gated cluster). Similar results are found for the case of gas 
expultion, as fenriched does not follow the predicted trend 
with GG mass. We also do not find any significant trend 
between fenriched and cluster metallicity, [Fe/H]. Our results 
are in agreement with the previous analyses performed by 
Garretta et al. (2010a) and Khalaj & Baumgardt (2015). 
We also note that GGs 1, 2, 3, and 5 in the Fornax Dwarf 
Galaxy have similar fenriched values (Larsen et al. 2014b). 

The constancy of fenriched as a function of Galactocen¬ 
tric distance and [Fe/H] shows that GGs belonging to the 
halo and the bulge do not differ in their present enriched 
fractions. This is relevant as metal poor halo clusters are of¬ 
ten thought to be accreted from dwarf galaxies while metal 
rich bulge clusters are thought to have formed in-situ (e.g., 
Brodie & Strader 2006). Hence, regardless of the formation 
epoch or location (and subsequent migration), GGs contain 
similar fractions of enriched stars. 

Our results suggest that GGs were not initially signifi¬ 
cantly more massive than at present, or at least did not lose 
their stars due to the influence of the tidal field or gas ex¬ 
pulsion (or birth location), contrary to the requirements of 
the self-enrichment scenarios discussed in § 1. This, in turn, 
implies that the “mass budget problem” cannot be solved 
by increasing the initial mass in clusters. 

If the clusters did not undergo a signihcant amount of 
preferential mass loss, then the fraction of enriched to pri¬ 
mordial stars observed in clusters likely reflects the initial 
value. If this is the case, then none of the self-enrichment sce¬ 
narios can work, without extreme IMF variations, it is not 
possible to produce the amount of material required to form 
the observed number of SG stars (even when considering the 
fraction of primordial material required to be mixed with the 
polluted ejecta) through standard nucleosynthetic channels. 
Additionally, if the current value of fenriched reflects the ini¬ 
tial value, the Galactic halo would not be made of significant 
numbers of lost FG GG star^. In turn, GGs, due to their 
lower initial masses, would have contributed substantially 
less to the reionisation of the Universe at high redshift than 
suggested in the heavy mass-loss paradigm (c.f., Schaerer & 
Gharbonnel 2011). 


^ Garretta et al. (2010a) find that 1.4% of the halo may be 
made up of enriched stars that have escaped GCs. A similar 
fraction has been found by Martell et al. (2011) and Ramirez 
et al. (2012). 
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4.2 Other Potential Solutions to the Mass Budget 
Problem 

Some recent models have attempted to solve the mass bud¬ 
get problem by having young GCs accrete enriched mate¬ 
rial from their surroundings, i.e., AGB processed material 
from other stars in the host galaxy not stars within the GC 
(e.g., Maxwell et al. 2014). However, this scenario is dis¬ 
favoured due to the lack of Fe spreads in most GCs observed 
to date, as the host galaxy ISM is expected to be enriched 
by both material processed through AGB stars, as well as 
high mass stars and their associated SNe products. Even 
small amounts of material processed in SNe should result in 
substantial Fe spreads, especially in low-metallicity clusters 
(e.g., Renzini 2013), so it is unlikely that the accreted mate¬ 
rial would not introduce (at least small) Fe spreads within 
the GCs. Additionally, the ISM will be polluted with the 
ejecta of lower-mass AGB stars, which do not match the 
abundance trends observed in GCs in general (D’Ercole et 
al. 2010 ). 

Alternatively, Charbonnel et al. (2014) have suggested 
a solution to the mass budget problem by forming the FG 
with a stellar IMF devoid of stars below 15 — 20 M 0 (i.e., the 
FG only consists of stars that contribute to the processed 
material needed to form the SG). Hence, all stars observed 
today would be SG stars, even the ones with “FG-like” abun¬ 
dances. Since there are no long-lived FG stars, having the 
FG be 10-100 times more massive than the current mass of 
the clusters would not contradict the current observations 
(unless dark remnants of the FG remained which would af¬ 
fect the M/L of the clusters, but these are assumed to be 
ejected due to SNe kicks). In this scenario, fenriched would 
represent the fraction of SG stars born with abundances ob¬ 
servably different than the SG stars born primarily from pri¬ 
mordial material. However, in this scenario large He spreads 
are expected (e.g., Chantereau et al. 2015), contrary to ob¬ 
servations (Milone 2015; Bastian et al. 2015) and no young 
massive clusters with such anomalous stellar IMFs have been 
observed to date (e.g., Bastian et al. 2010). 

We note that resorting to such an anomalously top- 
heavy stellar IMF is unlikely to help in the AGB scenario. 
Unlike in the FRMS scenario, the polluting stars do not 
end their lives as SNe, and hence will not receive strong 
kicks at the end of their lives that could potentially remove 
them from the cluster. Such remnants would leave a strong 
dynamical signature within clusters as well an anomalous 
population of white dwarfs, which would be inconsistent 
with observations of white dwarfs in some GCs (Richer et 
al. 2008). In the standard AGB scenario, the IMF is already 
assumed to be highly anomalous, as the SG only consists of 
stars that would remain alive until the present day, i.e. only 
stars with masses below 0.8 M©. Allowing for a normal IMF 
for the SG stars would result in an additional factor of two 
in the mass-budget problem (see Cabrera-Ziri et al. 2015 
for a further discussion of this point and other underlying 
assumptions made to minimise the mass-budget problem in 
the AGB scenario). 


4.3 Constraints from Other Sources 

Our results show that the observations are inconsistent with 
the proposed self-enrichment models, as GCs are unlikely to 


have lost substantial fractions of their initial masses. Similar 
conclusions have also been reached independently from a 
variety of sources. Along similar lines as the present work, 
Larsen et al. (2012; 2014a) have shown that GCs in a sample 
of three dwarf galaxies could not have been more than 5 — 
10 times more massive at birth than their present mass, 
in tension with the requirements of self-enrichment models. 
Similar results have also been found for the Galactic Bulge 
(Schiavon et al. in prep.). 

The above results also place constraints on the size of 
the initial GC population, namely whether it was much 
larger than observed today. Since the observed field stars 
in dwarf galaxies and the Galactic bulge are already not nu¬ 
merous enough to produce enough enriched material to form 
the observed enriched (second generation) stars, there is lit¬ 
tle room for a large initial GG population where the majority 
of GGs are completely disrupted. Hence the current GCs are 
not just the surviving minority (due to some special initial 
conditions) of an initially much larger GC population. 

The main reason why nuclear burning has been invoked 
to explain the origin of multiple populations, is that many 
of the elements observed or inferred to vary from star-to- 
star within GCs are associated with hot hydrogen burn¬ 
ing. However, in such processing, the abundance of He (the 
main product of hydrogen burning) is expected to be di¬ 
rectly linked with changes in other elements, such as Na, O, 
C, and N. Recently, Bastian et al. (2015) have shown that 
the observed extent of the Na-0 correlations in clusters is 
not directly linked to the He spreads within the clusters, 
in contradiction to basic nucleosynthesis. Additionally, they 
showed that if an enriching source can explain clusters like 
NGC 2808, it will necessarily fail to reproduce more typical 
GCs like NGC 104 (47 Tuc), and vice-versa, unless there 
is a strong stochastic element to the model. This is funda¬ 
mentally at odds with all the proposed enrichment sources 
(AGBs, FRMSs, very massive stars, or interacting binaries) 
and appears to rule out all nucleosynthetic sources as the 
origin of the abundance anomalies. 

All self-enrichment models proposed to date predict 
that the enriched stars should be more centrally concen¬ 
trated than the primordial stars (modulo dynamical evo¬ 
lution). However, the observed radial profiles of some GCs 
(e.g., M15 and NGC 2808) are inconsistent with predictions 
as the primordial stars are more centrally concentrated than 
the enriched stars (M15 - Larsen et al. 2015 ) or that all pop¬ 
ulations have the same distribution (NGC 2808 - lannicola 
et al. 2009; Dalessandro et al. 2011). These discrepancies 
remain even after considering the effects of dynamical evo¬ 
lution. 

Since the proposed scenarios do not invoke any spe¬ 
cial physics (i.e. conditions in the early universe) the pro¬ 
cesses should be observable in massive clusters forming to¬ 
day. However, to date, no ongoing secondary episodes of 
star-formation within these clusters has been found (Bas¬ 
tian et al. 2013a; Cabrera-Ziri et al. 2014) nor have the gas 
reservoirs required to form the secondary population been 
found (Bastian & Strader 2014; Cabrera-Ziri et al. 2015; 
Longmore 2015). 

Additionally, observations have found unexpected cor¬ 
relations between [N/Fe] and GC mass in the M31 GC pop¬ 
ulation (Schiavon et al. 2013) as well as a strong trend of He 
spread with GC mass in the Galaxy (Milone 2015). These 
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are not predicted in any of the self-enrichment models, and 
can only brought into agreement with models by making 
the mass-budget problem significantly worse (by additional 
factors of three or more - see Bastian et al. 2015). Finally, 
the discovery of the chemical anomalies in bulge stars of 
the Galaxy (that do not appear to be from dissolved clus¬ 
ters - Schiavon et al. in prep.) suggests that whatever pro¬ 
cess that causes the multiple population phenomenon is not 
limited to GCs, but rather may be a relic of high-redshift 
star-formation in dense environmnets. 

It is still unclear what the controlling factor is in de¬ 
termining whether a cluster hosts chemical anomalies. The 
standard idea is that cluster mass is the key parameter, al¬ 
though this was largely based on the idea that GCs were 
originally much more massive than they are today. The ob¬ 
servations presented here, along with studies of the stellar 
populations of dwarf galaxies (Larsen et al. 2012; 2014a) 
show that GCs are unlikely to have been signihcantly more 
massive in the past, hence their current masses (and fenriched) 
should approximately trace their initial values. Since young 
and intermediate age (< 2 — 3 Gyr) clusters with masses 
similar to that of GCs (< 2 x lO^M©) have not been found 
to display chemical anomalies (Mucciarelli et al. 2008; 2011; 
2014; Davies et al. 2009; Mackey et al. in prep.), it appears 
that GC mass is not the key parameter. 


5 CONCLUSIONS 

We have shown that the fraction of enriched stars in Galactic 
globular clusters is surprisingly uniform, fenriched = 0.68 ± 
0.07, with the scatter being nearly entirely accounted for 
within the error budget (both in numbers and the associated 
error of assigning if a star is enriched or not), leaving little 
room for correlations with other parameters. 

Scenarios that invoke heavy mass loss (preferentially 
from the FG) from young globular clusters in order to solve 
the mass-budget problem, have suggested that the mass-loss 
is due to tidal stripping, gas expulsion, or the influence of 
the birth environment. Such heavy mass loss is expected to 
result in correlations between the observed fraction of SG- 
to-FG stars and other cluster properties, such as their mass, 
metallicity, or Galactocentric distance. 

We have collected data from the literature for 33 GCs 
and no such correlations are found, suggesting that clusters 
have not undergone such radical mass-loss during their lives. 
The fraction of enriched stars is independent of the GC’s 
metallicity, Galactocentic distance, and mass, from which 
we infer that it is also independent of their birth location (in 
the Galaxy or accreted from dwarfs) or epoch of formation. 

Hence, the observed constancy of fenriched rnay re¬ 
flect the initial value, which would exclude self-enrichment 
through nucleosynthesis as a potential solution to the puzzle 
of multiple populations, as there would not be enough ’first 
generation’ stars to produce the material require to for the 
observed ’second generation’ stars. 

Given the results presented here, along with recent re¬ 
sults from the literature, it appears that the mechanism re¬ 
sponsible for the chemical anomalies observed in GCs is not 
self-enrichment, at least not through the proposed nucle- 
osynthetic channels. Therefore alternative theories need to 
be explored. 
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